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We investigated SU(3) lattice gauge theory with a fundamental and adjoint plaquette term in the action. The 
purpose is to test whether the choice of a negative adjoint coupling can reduce lattice artefacts and improve the 
scaling behaviour. To this end, we have studied the finite temperature phase transition, the static potential and the 
mass of the ++ glueball. We found that indeed the lattice artefacts in e.g. m ++ /T c can be reduced considerably 
compared with the pure Wilson (fundamental) gauge action at the same lattice spacing. 



1. INTRODUCTION 

The choice of the gauge SU(3) action plays an 
important role in QCD lattice simulations, and 
many efforts have been spent in searching for a 
formulation where discretisation errors are reduced 
and/or topological dislocations are suppressed. In 
this work we investigate a gauge action containing 
plaquettes in both fundamental and adjoint repre- 
sentation. The motivation for this choice is that 
in the (/3/, (3 a ) plane, the pure SU(3) gauge theory 
in four dimensions has a line of first order phase 
transitions with an endpoint at 



(/?/,&) = (400(7), 2.06(8)) 



(1) 



The presence of this transition line causes large 
lattice artefacts; in particular, the mass m Q ++ of 
the lightest glueball ++ is expected to go to zero 
as the end-point is approached |2] • As a relic of this 
behaviour, the estimate of m ++ in some physical 
unit at (3 a = 0, for 5.5 < f3f < 6.0 is much smaller 
than the continuum result. The purpose of this 
work is to study whether the discretisation errors 
can be reduced by choosing a negative value of (3 a , 
i.e. by moving away from the transition line. 
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2. THE ACTION 

We consider lattice SU(N) gauge theory with a 
mixed fundamental-adjoint action, 

i--ReTr/E/p + (2) 



S = f3 f 
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where (3 a are the fundamental and adjoint cou- 
plings and Up is the elementary plaquette. We 
adopted 4-dimensional lattices with spatial exten- 
sions aJV s and temporal extension aN t , with pe- 
riodic boundary conditions in all directions. The 
Yang-Mills theory is recovered in the naive con- 
tinuum limit if the bare coupling go is defined by 
6/ ff g = f3 f + 2f3 a . 

Details on the simulation algorithm and its per- 
formance are given in g|. 

3. THE FINITE TEMPERATURE 
PHASE TRANSITION 

For a lattice with N t points in the time direction, 
in the limit N s — > oo, the deconfinement tempera- 
ture is given by 

Y c =Ntai{fif,P*}c) : (3) 

where {Pf,f3 a } c indicates the critical coupling. In 
our study we kept fixed the adjoint coupling to 
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Table 1 

Numerical results for the finite temperature phase 
transition obtained for the fundamental- adjoint 
gauge action. 



N t ; (3 a 


0.0 


-2.0 


-4.0 


2 


5.0948(6) 


6.4475(6) 


7.8477(6) 


3 


5.5420(3) 


7.1603(3) 


8.8357(4) 


4 


5.6926(2) 


7.4433(3) 


9.2552(6) 


6 




7.8056(5) 


9.7748(11) 



(3 a = 0,-2.0,-4.0 and determined /3/, c for Nt = 
2, 3, 4, 6 by adopting the method discussed in g|. 

Our results for f3f tC are summarised in Tabled 
for f3 a = we find a substantial agreement with 
other results present in the literature |5l6l7l8j . For 
Pa = 0, N t = 6 we performed no own simula- 
tion, but used in the following the result f3f <c — 
5.89405(51) of ref. 0. 

4. THE STATIC POTENTIAL 

In order to fix the scale we computed the static 
potential (at zero temperature) and extracted the 
string tension a 2 a. The static potential has been 
computed through the Polyakov loop correlation 
function: 

aV(r)=~\\og{P(x)*P(y)) + e] , (4) 

where y — x + rl and e is the correction due to ex- 
cited states in the string spectrum. Here we used 
a large temporal extension aN t » 1/T C to ensure 
that finite- Nt corrections are negligible; in particu- 
lar we adopted N t = 6/ (aT c ) for all computations. 
The string tension has then been evaluated from 
the ansatz 

V{r)=or + v-JL(l + ±y (5) 

where b has been recently shown from theoretical 
principles to be zero [Hj- 2 Moreover, we adopted 
the tree-level improved distance n, defined such 
that the force at tree-level has no lattice artefacts. 
In order to compute the static potential up to large 
distances by keeping the statistical uncertainties 

2 In our evaluation we considered b = 0.04fm, which was 
evaluated numerically in 111)1 . Nevertheless the effect of 
having b ^ is only minor and not crucial for our compu- 
tation of a. 



under control, we adopted a variant of the algo- 
rithm proposed by Liischer and Weisz |1U|. The 
main difference consists in using factorisation in 
the spatial directions, in addition to the one in the 
temporal direction. The full details of the proce- 
dure are given in For l/(aT c ) = 4 we were able 
to extract a 2 a from the force at r/a = 7 (/3 a = 0); 
in the other cases our final values were taken from 
r/a — 4 and r/a — 6. In all computations we 
where confident the quoted error also covers pos- 
sible systematic uncertainties. The numerical re- 
sults for a 2 a are reported in Pj. 

The dimensionlcss quantity T c /^/a as function 
of the lattice spacing is plotted in fig. ^ together 
with other values obtained for the Wilson action 
[Sj. We notice that for f3 a = —2 and —4 the es- 
timate for T c /y/a is closer to the continuum limit 
than for (3 a = 0. However, the difference between 
X/(aT c ) = 3 and l/(aT c ) = 4 is larger than that 
for the different values of a at fixed l/(aT c ). 

5. THE 0++ GLUEBALL MASS 

The mass of the lightest glueball is expected to 
be most sensitive to the choice of the action in 
the (/3/,/3 a ) plane. The ++ glueball mass has 
been computed through the connected correlation 
function between spatial Wilson loops, by adopt- 
ing substantially the same method used in jll| . 
Also for the glueball 2-point function we made use 
of an error reduction procedure based on the idea 
proposed in [10) . and already applied for the com- 
putation of glueball masses ^2 ■ We extracted the 
glueball masses at t/a = 2, 3 for 1 / (aT c ) = 2,3 and 
t/a = 3,4 for l/(aT c ) = 4,6. Fig. Ushows our fi- 
nal results for m ++ /T c as function of (aT c ) 2 ; here 
the errors are dominated by the uncertainties on 
m ++. By averaging several results in the litera- 
ture |13I14I15I16| for m ++ro and then using the 
continuum limit relation jll| 

7> = 0.7498(50), (6) 

we obtain 

m 0++ /T c | a=0 = 5.73(9) (7) 

as estimation of the continuum result. At a ~ 
0.1 lfm we do observe a moderate reduction of the 
lattice artefacts by using j3 a < with respect to 
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Figure 1. Results for T c j \fa as a function of 
(aT c ) 2 . In addition, we report results from JJf for 
the Wilson action at smaller lattice spacings. 
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Figure 2. m ++/T c for j3 a = (filled squares) 
f3 a = —2 (open squares) and (3 a — —4 (triangles) 
as function of (aT c ) 2 . The circle gives the contin- 
uum results extracted from the literature. 



the usual Wilson action ({3 a — 0). For j3 a = 0, 
the deviation from the continuum result of eq. J7J) 
amounts to ~ 18%, while for (3 a = —2, —4 it 
slightly decreases to ~ 12%. 

At a ~ 0.1 7fm one observes discretisation errors of 
~ 40% for the Wilson action, while for the mixed 
action they amount to ~ 25% for j3 a = — 2 and 
~ 20% for /3 a = -4. 

6. CONCLUSIONS 

We investigated scaling properties of a SU(3) 
lattice gauge action with plaquette terms in the 
fundamental and in the adjoint representation, 
with negative adjoint coupling j3 a . By studying 
the scaling behaviour of the quantity T c /^/a for 
the different (3 a at our disposal, we did not observe 
a significant improvement at negative adjoint cou- 
plings in comparison to the Wilson case (3 a = 0. 
The values obtained with negative a are a lit- 
tle closer to the continuum limit. As expected, 
the mass m ++ of the lightest glueball is more 
sensitive to the variation of (3 a . We investigated 
the scaling behaviour of the dimensionless quan- 
tity m ++/T c . Here indeed, we observed a signifi- 
cant reduction of the lattice artefacts for negative 
j3 a . At a ~ 0.17 fm, the lattice artefacts for a = 



are 40%, while for (3 a — —4 they decrease to 20%. 
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